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ABSTRACT 

We present results of simulations of the spectrum of the accretion flow onto the su- 
permassive black hole in our Galactic Centre, Sagittarius A*, generated with a cou- 
pling of Monte-Carlo (MC) radiation and general relativistic magnetohydrodynamic 
(GRMHD) codes. In our modeling, we use the 2D HARM GRMHD code to first model 
the physical parameters of the disk, then feed its results into our 2D MC photon trans- 
port code. We will discuss results obtained which fit radio, IR, and Chandra-obtained 
flaring or quiescent x-ray data points, as well as the validity of the amount of scaling 
of input parameters (density, temperature, and magnetic field) required to fit these 
points. HARM output will be used to suggest whether the scaling is within reasonable 
limits. 
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1 INTRODUCTION 

In recent years, a great deal of effort has been made to un- 
derstand the complete picture of Sagittarius A* (Sgr A*), a 
radio/IR/x-ray emission source at our Galactic Centre. It is 
widely accepted that this source is related to the accretion 
flow of a supermassive black hole whose mass we have taken 



to be 3.6 million solar masses dGhez et al.||20Q3||Schodel et 
[gI][M7|peIia|[2006t . 

The spectrum of Sgr A* shows several important com- 



ponents or signatures. 'Balick & Brownj(1974) first discov- 
ered the source in the radio/NIR, and ensuing observations 
by a number of researchers confirmed the strength of the 
source to be primarily in these regimes. Years later, Sgr A*'s 
x-ray spectrum was explored by Baganoff et al. ( 2001| 2003 ), 
which found notable signatures higher than 10 Hz, and 
also found a variability in this range, which suggests sep- 
arate spectral states - flaring and quiescent. Simultaneous 
mult i- wavelength observations have generally supported this 
suggestion (Eck art et al.f2 004, 2006; Belan ger et al.f20 05). 

A number of models have arisen to explain the mech- 
anisms at work to produce the spectra of black holes. The 
standard-disk idea was first explored by [Shakura Sun-] 
yaev ([1973) and proposes a situation in which gravita- 
tional energy is efficiently converted to radiation. This is 
most appropriate for optically thick disks with the flows 
on nearly Keplerian orbits. Sgr A*'s luminosity is less than 
10~^ L^cZcZ, where the Eddington luminosity LEdd represents 
the luminosity level where the gravitational force inward 



equals the radiation force outward for spherical accretions. 
Sources with such low luminosities may be described by the 
standard-disk idea in cases where the disk remains cold and 
is optically thick, and the luminosity is low primarily due to 
a very low accretion rate. However, the accretion flow in Sgr 
A* must be very hot, and there is no evidence for a cold, op- 
tically thick disk component. The low luminosity is limited 
by both a low accretion rate and a low emission efficiency of 
the flow. In cases such as this, other accretion models need 
to be considered. 

One developed model is called the advection dominated 



accretion flow (ADAF) (Ichimaru 1977). This suggested that 
close to the horizon, much of the energy of the accretion flow 
was advected into the black hole, rather than being radiated 
away, and was explored by a number of groups, including 
Nar ayan Yi| ( [1994] ), [Abramowicz et ar |( |1995t , and [Yuan 
jet al.| ( |2003t . It was found that t he ADAF solution can fit 
Sgr A* data well ( Narayan et al.|1 995), but while promising, 
this approach had several major drawbacks, namely its one 
dimensional approach and simplification of magnetic fields. 
The X-ray emission is produced primarily at large radii, 
which cannot account for the observed short-timescale X- 
ray flares. In an attempt to consider an alternative approach. 



Yuan et al. (2003) suggested inclusion of nonthermal elec- 
trons and found satisfactory fits, but again with simple mag- 
netic fields. 

It became apparent that a more complicated treatment 
of magnetic fields would be important for a more accurate 
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simulation, as the magnetorotational instability was found 
to be of vital importance in the development of turbulence 
that drives the accretion flow fHawley Ba lbus|2002 Haw- 
ley||2001 Balbus 2003 ). This instability ensures that in an 
accretion disk environment the necessity of outward angu- 
lar momentum transport leads to the establishment of very 
complex flow and field patterns. 

Models have since begun using magnetohydrodynamic 
codes to simulate accretion flows, with some measure of suc- 



Ohsuga et al. (2005) ran emission simulations on an 



MHD model by Kato et al. (2004) and was able to fit data 
well, but had difficulties doing so without cutting out a large 
portion of their simulation volume. [Goldston et al.| ( |2005| 
calculated the radio spectrum based on MHD simulations 



done by [Igumenshchev et al. (2003), but not other spectral 
bands. 

We hope to show a simulation method which minimizes 
compromises on consistency via assumptions or simplifica- 
tions by using the best tools available to present a more 
realistic accretion disk picture. 

We will present our approach to Sgr A* modeling using 
GRMHD and MC methods in §2, present the results gained 
through these methods in §3, and revisit the most important 
revelations and conclusions made within in S4. 



2 SIMULATION METHOD 

Our modeling method involves the coupling of two very dif- 
ferent codes which complement each other to provide a con- 
sistent view of the physical and spectral conditions in the 
accretion disk of Sagittarius A*. Detailed in the appropriate 
sections below are descriptions of the codes, discussions of 
their appropriateness for our work, and the manner in which 
they were coupled to provide an overall model. 



2.1 HARM GRMHD physical space modeling 

The GRMHD code selected to determine the model's phys- 
ical space is the 2D HARM code presented and detailed in 



Gammie et al. ( 2003J and Noble et al. ( 2006 ) . It is not within 
the scope of this paper to completely detail the HARM 
code's inner workings, so readers should see these referenced 
papers for further information. 

The code solves hyperbolic partial differential equations 
in conservative form - uniquely suiting it to a number of 
astrophysical problems, specifically those involving magne- 
tohydrodynamics in areas where general relativity is impor- 
tant. As the code evolves the space through time, conserved 
variables are converted to primitive variables at each step, 
to calculate a set of fluxes, given a set of sources. Use of 
primitive variable calculations allows the code to work with 
analytic solutions, rather than finding solutions numerically 
- leading to much faster calculation time. 

In this case, the variables tracked include density, total 
energy, internal energy, flow velocities, and magnetic fields. 
The last two are calculated both as 3-component and 4- 
component tensor quantities. 



2.1.1 Model set-up for our work 

As imported, HARM was set-up quite appropriately for our 
purposes. It was configured to evolve an accretion disk about 
a black hole, given a number of user-controlled parameters, 
on a 2D grid spaced (in sperical coordinates) radially and 
angularly, and assumed to be axially symmetric about the 
black hole's spin axis. The cells are assumed to be uniformly 
spaced with regard to a set of coordinates XI and X2, which 
can be converted then to r and theta, respectively. This con- 
version leads to a logarithmic spacing in r, with cells closer 
to the horizon having greater resolution (and smaller size), 
and more concentrated cells closer to the equator. This ef- 
fectively increases resolution in the plane of the accretion 
disk and close to the black hole, where the detail is most 
useful due to the much smaller length scales of interest in 
these regions. The resolution increase toward the horizon is 
vital to the successful running of the code, as it helps main- 
tain proper cell aspect ratios on the polar grid. A graphical 
representation of HARM's grid layout can be seen in Figure 
5, where it is compared to the MC code's grid. 

As the start, HARM seeds an equilibrium torus around 
the black hole with density as detailed by [Fishbone Mon-| 
|crief| (|1976|). The torus is perturbed by adding a small 
poloidal magnetic field and allowed to evolve around the 
black hole. We have not attempted runs with an initial 
toroidal, azimuthal field, though this is allowed by axial sym- 
metry, and should serve to increase the total field at the end 
of simulation. The additional toroidal component should not 
influence the MRI development, instead only adding to the 
final field. 

While most default parameters as included are appro- 
priate for the simulation, several had to be tweaked for this 
project. To effectively simulate conditions near an a low lu- 
minosity Active Galactic Nuclei, where the nonrelativistic 
gas pressure presumably dominates, an adiabatic index of 
5/3 was chosen. As a first approximation, we have chosen a 
non-spinning black hole. Future work will study the effects 
of including black hole spin. The simulation volume ranges 
from just inside the horizon at 2 GM/c^ to 40 GM/c^ - with 
the initial torus having an inner radius of 6 GM/c^ and its 
pressure maximum at 14.7 GM/c^. Our final trial was done 
on a grid of 512x512 - that is, 512 radial cells by 512 angular 
cells. Output from the code is all scaled to M, the black hole 
mass, for near complete freedom. 



2.1.2 Results and interpretations 

Our simulation was run to approximately 8000 timesteps, 
which equates to about 1.4x10^ seconds, in physical time 
within the simulation. Figure 1 displays and explains four 
panels, which show density, internal energy, squared mag- 
netic field, and bulk lorentz factor through the simulation 
volume. Of importance to note is that the images are shown 
evenly spaced in cell number. This creates a kind of stretch- 
ing effect near the equator and near the horizon - causing 
the disk to appear much thicker than it actually is. 

As expected, the density concentrates itself around the 
equatorial plane, with the internal energy contours closely 
following. The magnetic field is strongest in the disk region, 
but saturates the region, generally being stronger nearer the 
horizon. More importantly, one can see how highly turbu- 
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Figure 1. Panels from final timestep t = 8000 (1.4x10^ s) of 
HARM data, as labelled for different physical quantities. We show 
a cut-out in the r-z plane, with the black hole's horizon shown 
on left. As detailed in the text, the radial direction is scaled log- 
arithmically, and angular dimension is scaled more finely toward 
the equator. Therefore, the finest resolution is obtained closest 
to the horizon and equator. This causes a "spreading" effect in 
the images, as they are shown evenly spaced by cell number, not 
physical value. Each quantity is shown on a logarithmic scale, 
with dark red being the highest values, black being the lowest. 

lent the field becomes after being seeded smoothly in the 
initial torus. As the turbulence is expected to provide the 
mechanisms for electron heating, in future work, the level to 
which the field is churned up will be important to analyze. 
The final panel shows the bulk gas lorentz factor, which is, 
as expected, closely related to distance from the black hole. 

Figure 2 shows profiles of density, temperature, and 
magnetic field along an equatorial slice of the simulation 
volume. Again of note is the high level of turbulence in the 
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Figure 2. Profiles of magnetic field, density, temperature, along 
the equator, as a function of radius. These quantities are normal- 
ized to the maximum value for each. 
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Figure 3. Profiles of optically thin synchrotron and 
bremsstrahlung emissivity. These values are normalized to 
the maximum for each set, so no comparisons should be made 
between the strength of the two components from this plot. 
These values are found by determining the emissivity within the 
closest cell and multiplying it by the radius and scale height at 
that radius. 

magnetic field. Within 5 GM/c^, the density and tempera- 
ture drop to very low levels. This corresponds well with the 
flow velocities as this is expected to be where the flow be- 
comes almost completely radial, and may be mostly a prod- 
uct of how quickly matter falls to the horizon at this point. 
Of some concern is how much these final profiles are due 
to the intial torus set-up. Further trials would be needed to 
compare initial torus location to final results. 

Figure 3 also shows a profile of the synchrotron and 
bremsstrahlung emission as a function of radius. To calculate 
this value, we multiplied the optically thin synchrotron emis- 
sivity and bremsstrahlung emissivity, respectively, by the ra- 
dius and density scale height at that radius. The synchrotron 
emissivity was found by averaging the magnetic field over a 
number of zones near the equator - this helps eliminate some 
of the high variability seen in the magnetic field. It should be 
noted that the overall shape did not change dramatically as 
this was averaged over a smaller or larger number of cells. It 
is apparent that the majority of emission emerges between 



4 Guy Hilhurn, Edison Liang, Siming Liu, and Hui Li 



Disk Thickness 




Figure 4. Profile of density along a constant radius r = 9 GM/c^ 
slice over theta, from pole to pole. For reference, at this radius, 
the scale height is 2.7 GM/c^. 
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for the observed variations in Sgr A*'s spectrum. As accre- 
tion variation is due not only to instabilities in the flow but 
also to changes in density of accreted matter, this factor of 
two is not considered a constraint on scaling. 



2.2 Spectrum determination with MC code 

The MC code used for photon emission and scattering has 
a long history and is discussed in a number of resources 
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dCanfield e _ t al. 1987, L iang & De rmer^l988 i Bottcher et al 
1998; Bott cher Liang 1998; Bo ttcher et ai.|2003| |Finke 
Bottcher,|2005| |Finke||2bQ7t ,Bottcher Liang||2001t . Fora 
complete treatment, readers should see these papers. 

In general, this code is a coupled MC/FP (Fokker- 
Planck) code. For our intents, the FP evolution of the elec- 
tron distribution was unnecessary at this stage, so it was 
turned off to allow a fixed temperature given by the HARM 
output. The code is set up on a 2D axially- symmetric cylin- 
drical grid, creating a (hollow or solid, depending on whether 
the inner radius is set to zero) cylindrical shape. Each zone 
is assigned a density, ion and electron temperatures, mag- 
netic field amplitude, and thermal and nonthermal distri- 
bution components. For our purposes, this is simply set to 
be a Maxwellian, but the code allows power law nonthermal 
distributions as well. The code allows emission from the vol- 
ume and boundaries, and emitted photons are tracked and 
allowed to scatter or absorb. 

This approach should be quite consistent, but it should 
be mentioned that the MC code does not include general 
relativistic effects. Close to the horizon, the photons' paths 
should be bent by the high gravity - but this code does not 
include this effect. However, as was shown in the previous 
section, densities, temperatures, and emissivities are quite 
low, relatively, near the horizon. Because of this, it is ex- 
pected that these GR curvature effects would be minimally 
important for computing the global spectrum. 



Figure 5. Mass and luminosity accretion rates onto the black 
hole through a surface at 28M as a function of HARM timestep 
- with each timestep being equal to approximately 17 physical 
seconds. 

radii of approxiately 8M and 28M. Areas inside and outside 
this region contribute decreasing amounts to the total emis- 
sivity. It is important to note the lack of contribution from 
the 28M-40M region to overall emission, as our simulation 
volume for the MC code only extends to 28M (as is discussed 
below, in relation to the overlap of our two simulation grids). 

To demonstrate actual disk thickness, density along a 
typical slice through a constant radius r = 9 GM/c^ is shown 
in Figure 4. 

Our last HARM output figure. Figure 5, shows the vari- 
ance of mass accretion onto the black hole through a surface 
at 28M through time. When comparing parameters required 
to generate spectra at quiescent and flaring data points, the 
amount that accretion varies could be important in deter- 
mining the appropriateness of the variation in values used 
between the two points. In general, at a late time, mass 
accretion varies up to a factor of two. The time scale of ac- 
cretion variation, on an order of several hours, is reasonable 



2.2.1 Coupling GRMHD output to MC input 

Shown in Figure 6 is an overlay of the grids for the two codes. 
It can be easily seen that regardless of how our MC grid is 
set up, it will almost invariably undersample and oversample 
the HARM grid at different locations. This is minimized as 
much as possible by making a very large MC grid - in our 
trials, we use a 95x95 vertical by radial grid. The grid can 
actually be quite a bit coarser than this before one notices 
a significant deviation in spectrum. As we cannot cover the 
entire region with a rectangular grid, we have chosen a box 
which is 28 GM/c^ radially and 56 GM/c^ vertically. This 
box is oriented so that its inner edge is at the horizon. 

As noted previously, the output from HARM is com- 
pletely scalable by black hole mass. For our trials, a rather 
well-accepted value of 3.6x10^ solar masses was assigned. 
We also have the freedom to choose maximum initial den- 
sity and have the other parameters scale with the density 
consistently. Density is scalable in this way since it only re- 
flects a change in accretion rate. 

In the HARM output, to determine temperature, we 
use internal energy. As almost all contributions to this come 
from ions, we can determine an ion temperature in each 
zone, but not an electron temperature. Coulomb coupling is 
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Figure 6. An exaggerated representation of the HARM and MC 
code grid overlays. The sphericahy spaced grid underneath is that 
of the HARM code. The rectangularly spaced grid overlaid is that 
of the MC code. Both grids are axially symmetric, with said axis 
on the left side in this image. 



very low due to the low densities expected in the accretion 
disk of Sagittarius A* (about 10^ to 10^° particles/cm^), 
and as the nonthermal heating and cooling mechanisms are 
not well understood, maximum electron temperature is kept 
as a free parameter for spectral fitting purposes in this pa- 
per. This means, in effect, we are setting a constant global 
ratio between T^, the temperature of the ions, and Te, the 
temperature of the electrons, and therefore assuming a two 
temperature flow, though this ratio is allowed to change be- 
tween trials. 

We consider the magnetic field given by the HARM data 
to be the MRI saturated field. As mentioned previously, had 
we input an initial toroidal field, this would have contributed 
to the final total field. For this reason, we use the MRI- 
saturated field as a lower limit to our field value and allow 
scaling above this. 

So, in effect, we set a maximum density, which in turn 
sets an ion temperature and lower limit to the magnetic field. 
The electron temperature is freely changed by changing the 
ratio between T^ and Te, and the magnetic field can be 
scaled up from its initial value as necessary. 



2.2.2 Trial and fitting procedure 

Figures 7-11 show observational data points for Sagittarius 
A*, whose distance is taken to be approximately 8 kpc, and 
trial fits. Circles shown represent data points, upside down 
triangles are upper limit values, and the bowties are flaring 
and quiescent points from Chandra observations presented 



arising from the acceleration of moving charges by a mag- 
netic field. When sufficiently energetic, these charges (elec- 
trons, in our case) produce a continuum spectrum whose 
flux and turnover frequency are directly related to the val- 
ues of electron temperature, magnetic field, and density in 
this region. 

In such a hot medium, synchrotron self-absorption be- 
comes important. This phenomenon occurs when a photon 
interacts with a charged particle in a magnetic field and 
transfers its energy to the particle. The low energy positive 
slope of the synchrotron curve in each of the trials shown 
below is due to synchrotron self-absorption. 

The difficulties arise in fitting the flaring or quiescent 
points in the x-ray simultaneously with the radio/IR data. 
The approach made in this paper is to initially try to fit the 
flaring point. Its flat slope in L^; suggests that it might be 
fit by several possible components of the spectrum: 

(i) Bremsstrahlung emission, or free-free emission, arises 
(usually, and in this case) from the acceleration of a free 
electron by a free nucleus, in a completely ionized plasma. 
This creates a relatively flat spectrum in out to a cut- 
off point which corresponds with the temperature of the 
plasma. The flux of this component is directly related to 
the square of the density, as it depends on the population of 
electrons and ions, and also to the square root of the tem- 
perature. Temperature also serves to locate the high energy 
cut-off in the spectrum. Magnetic field values do not affect 
bremsstrahlung emission. 

(ii) Compton scattering is produced when photons inter- 
act with particles (in this case, electrons), leading to changes 
in energy for both photon and particle. Of interest for our 
work is inverse Compton scattering, where a less energetic 
photon gains energy in an interaction with a hot electron. 
In general, this tends to form a photon population whose 
shape is related to the initial photon distribution and the 
electron distribution it scatters from. This spectrum is typi- 
cally shifted up in energy an amount approximately equal to 
the square of the Lorentz factor of the electron population, 
and drops in flux by an amount equal to the optical depth 
of the scattering medium. For our trials, the Compton spec- 
trum of importance is mostly generated via synchrotron self- 
Comptonization (SSC), which refers to the Comptonization 
of a photon spectrum produced via synchrotron emission 
by the population of electrons responsible for emitting it to 
begin with. As the final spectrum resembles the electron dis- 
tribution and initial (unscattered) synchrotron spectrum, it 
is possible to locate this spectrum so its flat top intersects 
the flat flaring x-ray point. 

(iii) The second Compton scattering component, as its 
name suggests, is the spectral component created by the 
Compton scattering of an already scattered photon popu- 
lation. In this case, it shares the characteristics of the first 
scattered spectrum, and therefore is a viable component to 
fit the flaring x-ray point. 

All information on emission processes can be found in 



Rybicki & Lightman (1979). 



by Baganoff et al. (2001). 



Data in the radio to IR range are fairly easy to fit on 
their own. This is typically done with synchrotron emission 



A large number of MC runs were done to find several 
sets of parameters which fit this point while remaining con- 
sistent with the order of magnitu de estimate f or total lumi- 
nosity of 10"^^ erg/s suggested by 



Yuan 



(2007). 



The second step in fitting was to attempt to fit the qui- 
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Figure 7. Fit to flaring data with bremsstrahlung component. 
Maximum values of the 95x95 cell grid are scaled to n = 4.3x10^ 
particles/cm^, T = 1.2x10^ keV, B = 9.9x10^ G. Total luminosity 
is 3.4xl036 erg/s. 



Figure 8. Fit to flaring data with second Compton bump com- 
ponent. Maximum values of the 95x95 cell grid are scaled to n = 
3.4x10^ particles/cm3, T = 8.2x10^ keV, B = l.GxlO^ G. Total 
luminosity is 7.5x10^^ erg/s. 



escent x-ray data point by varying parameters from those 
found to fit the flaring point. This was done, specifically, 
by using the same input values as the flaring point, but by 
dropping either the density or temperature until the spec- 
trum went through the quiescent point. Depending on the 
component being fit to the x-ray data, one or both of these 
parameters could be changed to lower the luminosity to in- 
tersect the quiescent point. As we have data showing the 
variability of mass accretion by HARM, we can estimate 
reasonable changes in density. Due to the unknown nature 
of electron heating, it is not inconsistent to suggest that 
flaring and quiescent points may also arise due to brief mo- 
ments of higher or lower heating, caused by dissipation of 
large current sheets in the MRI-turbulence cascade. 

It should be noted that observations show a strong cor- 
relation between x-ray and radio/IR flux during flares. As it 
is generally accepted that the lower energy spectral com- 
ponent is due to synchrotron, this lends strength to the 
idea that the x-ray data is of synchrotron origin - either 
as a self- Comptonizat ion component, or an extended syn- 
chrotron component (|Liu &: Meiia|2QQl(|Markoff et al.|2Q01 



Dodds-Eden et al.|2009 ) 



3 RESULTS 

The results presented are shown first for fits to the flaring 
data (the higher x-ray point), then for fits to the quiescent 
data (the lower x-ray point). The HARM data used in the 
MC runs are based on output at the last timestep (t=8000). 
Fits are described by the component used to fit the x-ray 
data. However, this does not indicate that other spectral 
components were turned off during these runs. For instance, 
the bremsstrahlung spectrum can be seen to have both first 
and second Compton bumps, but the spectrum at the x-ray 
point is almost purely bremsstrahlung. 



3.1 Flaring results 

Figure 7 shows an attempt to fit the flaring data with the 
spectrum's bremsstrahlung component. This is similar to 



have found that the bremsstralung component is too high 
when the radio data are fit well by the synchrotron data. 
In this case, we h ave a simulatio n volume which extends to 
28 GM/c^ while lOhsuga et al.l (^2005) has volumes which 



extend to either 10 Rs or 30 Rs (20 GM/c^ or 60 GM/c^ 
respectively) . They found adequate fits with the smaller vol- 
ume with a density maximum of 1x10^° particles/cm^ and 
with the larger volume at 3.4x10^ particles/cm^. As would 
be expected, our trial lies somewhere between the two. How- 
ever, it seems unlikely that any skewing of density or tem- 
perature will allow a fit to the slope of the quiescent data. 

Figure 8 shows a fit by lowering density from the 
bremsstrahlung fit, but compensating by increasing temper- 
ature until the flat portion of the twice-scattered Compton 
bump fits the slope at the flaring point. Total luminosity 
and the spectral slope at the flaring point suggest this is a 
good fit - and the shape of the spectrum seems to allude to a 
possible fit to the quiescent point's slope if density is allowed 
to drop. As with the bremsstrahlung fit, we can compare the 
density maximum used in this trial to those of jOhsuga et| 



al, 



fits attempted by Ohsuga et al. (2005|. Like their trials, we 



2005| . The value here correlates to that group's lower 
density and higher volume value - this is reasonable, as we 
have a significantly higher temperature now. 

Figure 9 shows a fit with low density and high temper- 
ature, so the first Compton bump fits the fiaring x-ray data 
point. The synchrotron/synchrotron self-Compton approach 
was used analytically by |Liu et al.| ( |2006| ), which, while a one 
zone approximation, found very similar values to our max- 
imums for one of their trials: n = 2.3x10^ particles/cm^, T 
= 1.1x10^ keV, B = 10.2 G. This fit, like the one above, is 
promising due to its shape. By lowering only temperature 
from this value, it appears the quiescent point may be fit 
with a similar slope. 



3.2 Quiescent results 

As alluded to above, quiescent fits for the bremsstrahlung 
trial, second Compton trial by varying the temperature, and 
first Compton trial by varying the density could be found 
which intersect the point. However, all of these had nearly 
zero slopes in L^;, and were therefore unsatisfactory. We 
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Figure 9. Fit to flaring data with flrst Compton bump compo- 
nent. Maximum values of the 95x95 cell grid are scaled to n = 
5.2x10^ particles/cm^, T = 1.8x10^ keV, B = 13 G. Total lumi- 
nosity is 2.6x10^^ erg/s. 



Figure 11. Fit to quiescent data with flrst Compton bump com- 
ponent. Maximum values of the 95x95 cell grid are scaled to n 
= 5.2x10^ particles/cm^, T = 4.7x10^ keV, B = 13 G. Total 
luminosity is 3.4x10^^ erg/s. 
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Figure 10. Fit to quiescent data with second Compton bump 
component. Maximum values of the 95x95 cell grid are scaled to 
n = 6.9x10'^ particles/cm3, T = 8.2x10^ keV, B = 1.6x10^ G. 
Total luminosity is 8.2x10^^ erg/s. 



present the two quiescent fits which appeared most promis- 
ing. 

Figure 10 shows a quiescent fit using the second Comp- 
ton trial shown in Figure 8. It should be noted, first of all, 
that the choppiness in this spectrum is due to lack of statis- 
tics in the twice-scattered photons. Optical depth directly 
along the equator in this case is approximately lxlO~^, lead- 
ing to a very small population of scattered photons. The 
slope at the quiescent point is significantly steeper than the 
previous trial at the flaring point, and fits well the slope 
required at quiescence. In this case, we have dropped the 
density by less than a factor of five. This is slightly higher 
than the accretion variation of HARM data, a factor of two, 
but is not extreme. As noted in the variation analysis, vari- 
ations in the density of accreted matter are not considered 
by HARM, and therefore, changes of this level are not pre- 
cluded. The results shown here can be changed greatly by 
small changes in input values, due to the combination of 
spectral components required to produce this spectrum. 

Figure 11 shows the quiescent fit for the first Compton 
trial shown in Figure 8. This fit was obtained by dropping 
the temperature by about half an order of magnitude. As 



noted before, without better understanding of acceleration 
mechanisms, it is difficult to say whether global temperature 
changes of this scale are feasible or not. However, by this 
trial, this approach looks promising, and is closely in line 
with the idea that correlated low and high energy flares arise 
due to synchrotron radiation - in this case, the x-ray being 
fit by an SSC component. This trial gives our best fit to the 
quiescent data point's slope. 



4 CONCLUSIONS 

We have presented results for coupled physical and spectral 
simulations of the accretion disk of Sagittarius A*. Fits to 
flaring x-ray data are shown for three sets of parameters of 
an entirely thermal plasma, and these results are discussed, 
noting their similarities to two other recent models made 
for this same situation. The values used in our trials for all 
three sets of parameters are quite similar (respectively) to 
those used by L iu et al. (2006 ) and Ohsuga et al. ( |2005| ) . 

Of these, the trials which fit the flaring data point with 
the first Compton bump and second Compton bump show 
promise to allow fits to the quiescent x-ray data by only 
changing one physical value. The second Compton bump 
trial shows a promising fit which intersects the quiescent 
point by dropping the density, and may be at an appropri- 
ate slope. This trial is also the best fit to low energy data 
while fit to the x-ray points. While intersecting each x-ray 
point, some of the IR data are fit by the synchrotron bump. 
The first Compton bump trial shows a very good fit, both 
in slope and position, to the quiescent point by dropping 
the temperature, though these trials never match the IR or 
radio data points. The second Compton bump fit appears 
approximately valid in comparison to the variations in mass 
accretion onto the black hole suggested by the HARM out- 



put. It should be noted, as reported in Xu et al. (2006), the 
quiescent state X-ray emission can be attributed to thermal 
emission from the large scale accretion flow. For this reason, 
the data point we have fit here could be considered an upper 
limit. 

Overall, we present our second Compton bump trials as 
our best explanation for the observed spectra of Sagittarius 
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A*. Due to the combination of spectral components at the 
critical point in the x-ray, we are able to describe a large 
number of different slopes and locations - with these being 
quite sensitive to input values. This may account for the 
observed variability in x-ray slopes and luminosities. 

While our second Compton trial can fit the IR points, it 
is found that none of our fits can match the radio data well 
while also fitting x-ray data - this suggests that a completely 
thermal approach using GRMHD and one-temperature elec- 
trons is lacking in ways that may be improved by using FP 
techniques or results suggested by recent kinetic simulations 
of turbulence- heated electron distributions (lLiangll2009l). It 
is also possible that a larger simulation volume must be en- 
closed to fit these points concurrently with the x-ray data. 

In the future, this project will turn to some method of 
consistently introducing nonthermal particle acceleration to 
determine whether this can provide better fits to the joint 
radio/IR/x-ray data. 



ACKNOWLEDGMENTS 

This work was supported by a Los Alamos National Labo- 
ratory IGPP research grant and NSF-AST-0406882. 

The authors would like to thank a number of people 
whose input helped shape this work, including, but not lim- 
ited to, Markus Bottcher, Xuhui Chen, Justin Finke, Charles 
Gammie, and Feng Yuan. 



Hawley J. F., Balbus S. A. 2002, ApJ, 573, 738 
Ichimaru S. 1977, ApJ, 214, 840 

Igumenshchev I. V., Narayan R., Abramowicz M. A. 2003, 

ApJ, 592, 1042 
Kato Y., Mineshige S., Shibata K. 2004, ApJ, 605, 307 
Liang E., 2009, arXiv: 0902.4740 
Liang E. P., Dermer C. D. 1988, ApJ, 325, L39. 
Liu S., Melia F. 2001, ApJ, 561, L77 

Liu S., Petrosian V., Melia F., Fryer C. L. 2006, ApJ, 648, 
1020. 

Markoff S., Falcke H., Yuan F., Biermann P. L. 2001, 379, 
L13 

Martin K. W., Liu S., Fragile C, Yu C, Fryer C. L. 2009, 

arXiv: 0904.0118 
Melia F. 2006, The Galactic Supermassive Black Hole 

(Princeton: Princeton Univ. Press) 
Narayan R., Yi I. 1994, ApJ, 428, L13 
Narayan R., Yi I., Mahadevan R. 1995, Nature, 374, 623 
Noble S. C, Gammie C. F., McKinney J. C, Del Zanna L. 

2006, ApJ, 641, 626. 
Ohsuga K., Kato Y., Mineshige S. 2005, ApJ, 627, 782. 
Rybicki G. B., Lightman A. P. 1979, Radiative Processes 

in Astrophysics (New York: Wiley) 
Schodel R., et al. 2002, Nature, 419, 694. 
Shakura N. I., Sunyaev R. A. 1973, A&A, 24, 337 
Xu Y., et al. 2006, ApJ, 640, 319 

Yuan F., Quataert E., Narayan R. 2003, ApJ, 598, 301 
Yuan F. 2007, ASPC, 373, 95. 



REFERENCES 

Abramowicz M. A., Chen X., Kato S., Lasota J. -P., Regev 

O. 1995, ApJ, 438, L37 
Baganoff F. K., et al. 2001, Nature, 413, 45 
Baganoff F. K., et al. 2003, ApJ, 591, 891 
Balbus S. A. 2003, ARA&A, 41, 555 
Balick B., Brown R. 1974, ApJ, 194, 265 
Belanger G., Goldwurm A., Melia F., Ferrando P., Grosso 

N., Porquet N., Warwick R., Yusef-Zadeh F. 2005, ApJ, 

635, 1095 

Bottcher M., Jackson D. R., Liang E. P., 2003, ApJ, 586, 
339 

Bottcher M., Liang E. P. 1998, ApJ, 506, 281. 
Bottcher M., Liang E. P. 2001, ApJ, 552, 248. 
Bottcher M., Liang E. P., Smith I. A. 1998, A&A, 339, 87. 
Canfield E., Howard W. M., Liang E. P. 1987, ApJ, 323, 
565. 

Dodds-Eden et al. 2009, arXiv: 0903.3416 
Eckart A., et al. 2004, A&A, 427, 1 
Eckart A., et al. 2006, A&A, 450, 535 

Finke J., August 2007, Physics and Astronomy, Monte 
Carlo/Fokker- Planck simulations of Accretion Phenom- 
ena and Optical Spectra of BL Lacertae Objects 

Finke J. D., Bottcher M., 2005, PASP, 117, 483 

Fishbone L. G., Moncrief V. 1976, ApJ, 207, 962. 

Gammie C. F., McKinney J. C, Toth G. 2003, ApJ, 589, 
444 

Chez A. M., et al. 2003, ApJ, 586, L127 
Goldston J. E., Quataert E., Igumenshchev I. V. 2005, ApJ, 
621, 785 

Hawley J. F. 2001, ApJ, 554, 534 



